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It is well known that a spin-1/2 antiferromagnetic chain with nearest-neighbor coupling 
Ji and next neighbor coupling J 2 could dimerize in the presence of frustration, as exemplified 
by the exactly soluble one- dimensional Majumdar- Ghosh (MG) model ][]], where the exact 
(twofold degenerate) ground state is a simple product of singlet dimers. The elementary 
excitation can be constructed as a pair of unbound spins above the completely dimerized 
state 0. This model also can be considered as a two-chain with one diagonal or a zigzag 
spin ladder. Analytical and numerical studies of the MG model (or models alike) || show 
a transition from a gapless phase when J2 < J2C to a gapped phase when J2 > J2 C - In 
addition, a family of spin-ladder models were studied by Kolezhuk and Mikeska |3J] that 
exhibit non-Haldane spin-liquid properties as predicted by Nersesyan and Tsvelik H . Great 
progress were made in the studies of the ladder systems 0. In two dimensions, Shastry 
and Sutherland proposed a model (SS) with dimerized eigenstate twenty years ago [0]. The 
model had its experimental realization in compound SrCu2 (1^03)2 recently ||. Theoretical 
and experimental investigations of the SS model are currently undergoing ||. 

In this Letter we consider a set of quantum spin models, we call it the net spin model, 
with any spin S defined on a double layer as shown in Fig. 1. Each layer has M x L sites 
and it connects to the other not only perpendicularly by coupling J 1( but also diagonally by 
J3 and J4. L is a measure of the dimensionality varying from 1 (ID, the net spin ladder) to 
M (2D, the net spin layer). Throughout this work, we use S to represent spins on the lower 
layer, and S' on the upper layer. Open boundary conditions are imposed. 

Under certain combination of coupling constants, we show that for general L and M 
and any spin S the model could be solved exactly for a subset of its Hilbert space. First, 
depending on parameters, the model have two kinds of ground state, one is the completely 
dimerized state, and the other one is the ground state of the spin-a (a = 2S) model defined 
on the single layer. Second, the model exhibits rich excitation phases, depending on coupling 
constants and dimensionality. One may think of the net spin model as a generalization of 
the Majumdar- Ghosh model for any spin S in two dimensions, with much more rich phase 
diagrams than that of the MG model. 
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The net spin model Hamiltonian is: 



M,L 

H = 2Ji Sk t i ■ S'kj 

k,l=l 
M-l,L-l 

+ 2J 2 X! ^k,i • (Sfc+i,/ + Sfc,z+i) + S'k,i ■ (S'k+1,1 + S'k,i+i] 

k,l=l 
A/-1,L-1 



+ 2J3 S'k,i ■ (Sjfc+i,/ + Sfc^+i) 

k,l=l 
M-X,L-l 

+ 2J 4 Sk,i ■ (S'k+1,1 + S'fc^+i) , (1) 

k,l=l 

with Ji, J 2 , J3, J4 > 0, and |S| = |S'| = S. Both Jj and J2 favor local antiferromagnetic 
ordering while J3 and J4 represent frustration effects. 

For any integers M and L, we found that a complete dimerized state: 

Vfc = [1, 1'] % 2'] • • • [M, M'] ■ ■ • [N, N'], (N = ML) (2) 

is an eigenstate of Eq. (|T|) when 

2 J 2 = J 3 + J 4 • (3) 

To prove our statement, we first show that this is true for any 4-spin plaquette, 
{S^i, S'i,!, S2,i, S'2,i}, simplified as {Si, S'i, S2, S' 2 }, as indicated in Fig. 1. Since ipo is 
a direct product of singlet [1, 1'] and [2, 2'], we have 

(Si + S\)iP D = (S 2 + S' 2 )^d = , (4) 

H = J 1 [(S 1 + S'1) 2 + (S 2 + S' 2 ) 2 - (S 2 + S' 2 + S 2 + S' 2 ) " 

+ 2J 2 (Si + S'i) • (S 2 + S' 2 ) + (J 3 + ^ - 2J 2 )(S 2 ■ S\ + Si • S' 2 ) . (5) 

Immediately we see that dimerized state t/>d is an eigenstate of H with eigenvalue —AJ 1 S(S+ 
1) when condition Eq. @) holds. 

In Eq. (H), [i,j] denotes the normalized wave function of spin singlet |Sj + = 0, with 
Clebsch- Gordon coefficients for spin-S 1 dimer given by: 



We know that the formation of dimers [1, 1'] and [2, 2'] is due to antiferromagnetic cou- 
pling J\. Without antiferromagnetic couplings J 3 and J 4 , quantum fluctuation due to J 2 will 
kill the formation of both dimers. As long as condition Eq. @ holds, the effect of J 2 on the 
dimers will be cancelled out exactly by J3 and J4. Thus, no matter how many plaquettes 
we put together to form the two dimensional lattice as shown in Fig. 1, ip£> is the eigenstate 
of Hamiltonian Eq. (|IJ) with eigenvalue E D = -2J i iV5'(S' + 1), N = LM. 

It is also obvious that the dimerized state ipD is the ground state in the limit of J\ — > 00, 
and in the limit of J\ — > 0, ipo is not the ground state, so there exists a critical value J\ c 
such that when J x > J lc the ground state is completely dimerized. 

We should mention that for the case of L = 1, S = 1/2, the net spin model (the net spin 
ladder) goes back to the MG model |]] when J 2 = 0.5Ji, J 3 = Ji, and J 4 = 0. So our model 
is a generalization of the MG model for any spin S . Moreover, our model applies to two 
dimensional geometry and its ground state and the first excited state are exactly soluble for 
a wide range of parameters. 

To study the quantum phase transition and the excitation spectrum of the net spin model, 
we concentrate on two cases, one is L — 1, corresponding to the one- dimensional lattice, 
and the other one is L = M, corresponding to the two-dimensional lattice. We present 
qualitative discussions for cases 1 < L < M at the end. For the purpose of illustration, we 
take J 3 = J 4 = J 2 so that it enables us to get analytical solutions. Such simplification does 
not alter general features of the model, although to get exact numbers, one needs to perform 
numerical calculations for other combination of parameters { J 2 , J3, J4}. For L = 1, several 
studies on the generalization of the MG model, P,|3|JT0|, were carried out recently. 

We start from the ladder case, L — 1, the Hamiltonian can be rewritten as 

M M-l 

H = -2MJ X S(S + 1) + Jx £ *\ + 2J 2 °k ■ <Tk+i , (7) 

fc=l k=l 

where <r fc = + SV This Hamiltonian describes a chain of M spins with |er fc | range from 
(singlet, dimer), 1 (triplet), • • •, to the highest angular momentum 25'. 



With this Hamiltonian, it is easy to see that ipo is an eigenstate of H by using the 
relation (S& + S'k)[k, k'] = 0, hence 

/ M M-l \ 

U Y, °i + 2 ^2 E *k ■ >= , (8) 

V k=l k=l ) 

and the eigenvalue is E D = —2MJiS(S + 1). 

To lower the energy, the squared term, J\cr 2 k1 favors the system to be in the |cr fc | = dimer 
state, while the exchange term, J2(Tk • o~k+i, favors the system to be in the higher |tr&| state. 
If all pairs + S'& are in the same angular momentum a state, then the model is equivalent 
to the spin-cr chain with coupling constant 2J 2 , apart from a constant proportional to MJ\. 
One can see from Eq. ([7|) that the eigenstates of the spin-cr chain is also the eigenstate of 
the net spin ladder model. In fact, as J\ decreases further, the ground state of the net spin 
ladder model becomes the ground state of the spin-cr chain. The transition point (J 2 / ' J\) c 
is determined by 

m / J \ M_1 

(E^ + 2 t E^-^ + i) = o, (9) 

fc=l Wl/ c k=l 

where the expectation value is taken in the ground state of the spin-cx chain. To evaluate the 
ground state energy of the spin-a model, we observed that the linear spin wave theory [TT 



gives the ground state energy of the spin chain within 3% accuracy as compared with exact 
calculations such as the exact diagonalization |12|] and the density matrix renormalization 



group method [fL3|1 . So even though it is inappropriate to describe the magnetic properties 
such as the staggered magnetization in one and two dimensions, it is perfectly safe to use 
the spin wave theory for the estimation of the transition points (J 2 / 'Ji) c , given by 

(-) = ~ 7 a + ! " (10) 
\jj c 2vr(cr + l)-2 V ; 

For spin S — 1/2 and 1, we have verified those results by exact diagonalization calculations. 

Furthermore, we can show that once (J2/J1) > (J2/ Ji)c, the ground state jumps from 
all crfc = (the dimer state) to all = 2S, without passing any intermediate values of a^. 
The transition is of the first order. 



We next study the excitation spectrum of the model. Let us establish an important 
fact first. Instead of using {Sk,S'k,k = 1,2, ■ ■ ■ , N = LM} as basis operators, one can 
also use {cr fc = Sfc + S'k,k = 1,2,---,N = LM} to specify eigenstates. The complete 
dimerization state would be written as ipn = [0][0] ■ ■ • [0] . Then, one can show that ip m = 
[0] • • ■ [0] [m] [0] ■ ■ • [0] , and ijj m ,m' = [0] • • • [0] [m] [0] [m'] ■ ■ ■ [0] , and wave functions consisting 
nonzero as separated by at least one dimer (<Jk = 0), are also eigenstates of the Hamiltonian 
Eq. (§. 

For sufficient large value of J\, the ground state is the complete dimerized state and one 
can create a single magnon simply by changing one from to 1. For example, for spin- 1/2 

case, 

^ m = [l,l'][2,2']---(T,T)---[M,M'], (11) 

is the eigenstate of the net spin model with total spin S to t = 1- The energy gap in this case 
is exactly A st = 2Ji and the transition is from singlet to triplet. This excitation energy is 
independent of other antiferromagnetic couplings as long as they satisfy the condition Eq. 
(0) and thus it is possible that there exist states with lower excitation energy. 

Indeed, as we have shown in detail for the spin-1/2 case [|14|| , when (J2/J1) > 1/2, there 
exists another singlet with energy lower than the triplet. The corresponding wave function, 
defined as ip n with eigenvalue E s = Ji — 4J 2 , is that all pairs are dimerized except on 
one plaquette where the two dimers are broken into two triplets to make up a non-dimer 
singlet. This state is also highly degenerate for this plaquette could be anywhere. We can 
show that when 2J 2 > Ji, two isolated triplets embedded in otherwise all dimers would have 
higher energy than that of combining the two triplets next to each other to form a singlet. In 
another words, there exists an attractive interaction between the two triplets. Consequently, 
as more and more triplet excitations are created, a kind of phase separation would occur 
where the system consists of clusters of triplet and singlet. Note that this critical value 
(J2/J1) = 1/2 is independent of lattice size and dimensionality. 

For the general spin-S 1 case, same conclusions hold. Using the fact that ip m (E t ), ip n {E s ), 
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i>m,m' (2-Kt), etc., are also eigenvectors, we can show that 

A st = E t - E D = 2 Ji (12) 
A ss = E S -E D = 4( Ji - J 2 ) (13) 
A ss < A rt , J 2 /Ji > 1/2 . (14) 

To be more specific, we use S = 1 as an example and summarize our findings in Figure 2a. 
There are three regions, depending on parameter J 2 /Jj_: 

Region I , J%j J\ < 1/2, A/Ji = A st /Ji = 2, independent of J 2 and it is singlet to triplet. 

Region II , 1/2 < J 2 /Ji < (^2/^1)0 ^/J\ = A^/Jx = 4(1 — J 2 /Ji) decreases linearly 
from 2 to 4(1 — ( J 2 / Ji)c) = 1.4760 at ( J 2 / J\) c = 0.6310. The excitation is from singlet 
(dimers) to another singlet. In both regions I and II, the completely dimerized state 
ipD is the ground state. 

Region III , J 2 / J\ > ( J 2 / Ji) c , the ground state is no longer the completely dimerized state 
ipD- Rather, it is the ground state of spin a = 2 chain with coupling 2J 2 , degenerate 
with ipD right at (J 2 /Ji) c . The excitation is the Haldane gap of the spin-2 chain, 
A/ Ji ~ 0.088 [0 If the original spin is S, then in this region the ground state is that 
of the spin-2S" chain. For large S, the Haldane gap is approximately A pa 9.5S 2 e~ 7T,s . 

We now turn to the double layer case (L = M and iV = LM). Similarly, we get 

N N 

H = -2^5(5 + 1) + Ji "I + 2 J 2 E *k ■ vi, (15) 

k=l <k,l> 

where <Tk = Sf-j + S'k,i and < k,l > refers to nearest neighbors in two dimensions. Again, 
when all |cr fc | are the same, this is nothing but the two-dimensional spin-a antiferromagnetic 
Heisenberg model with coupling constant 2J 2 , apart from a constant. The critical point at 
which the complete dimerized state becomes the excited state is determined by, 
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and the spin wave theory gives 

J 2 \ 1 a 



(17) 



JJ C 2 2(7 + 0.3159 

As one increases J2, similar to the chain case, the change of Ck in the ground state is directly 
from (Tfc = to cr/c = 2S, no intermediate values of appear. 

For the excitation spectrum, again we can show that it is still true when J2/J1 < 1/2, 
the first excitation is a triplet. However, for the layer case, without detailed calculation of 
the ground state energy of the corresponding spin-cr AFH model, one can easily show that 
{J2I Ji)c < 1/2- This is because on the 2D square lattice, due to quantum fluctuations, the 
true ground state energy is always lower than that of the Neel state, —2Na 2 . Moreover, 
when (J2/J1) > {J2/ J\)ci the system is gapless because the existence of long-range-order in 
the 2D square lattice spin-cr Heisenberg model for a > 1, as shown rigorously long time ago 
Il^ , |T7j| . Therefore, there exists no transition from the dimer singlet to another singlet as 



what we have seen for the net spin ladder. We plot the excitation gap as function of J2/J1 
in Figure 2b. 

What will happen for cases of 1 < L < 00? Exact estimation of the ground state energy 
of spin-cr model is not necessary. As we have shown, whether the excitation spectrum will 
be ID or 2D like depends on whether (J2/ 'Ji) c > or < 1/2. According to the linear spin 
wave theory, we can write 

<!>».,,> = i^±. (18) 



<k,l> 



where a = 2 — 1/L accounts for the Neel state energy, and b ranges from 0.3634 (ID) to 
0.3159 (2D). Exact account of the quantum fluctuation will not change this number much. 
We can safely use these numbers to make estimations. For spin S = 1/2, we have studied 
L = 2 and L = 3 cases [TJJ] by using the spin wave theory, perturbation theory, and exact 
diagonalization on small lattices, they are all consistent. Results indicate that there seems 
to have a dimensional crossover from ID to 2D occurred at L = 3 to L = 4. 

For other combinations of couplings, one can only determine the critical points and 



eigenvalues of the model numerically (see, for example Ref. |T8|). However, qualitative 



behaviors should be the same as what we have shown here. Note that the energy gap 2 J\ is 
quite large (critical values of J 2 are smaller than Ji) so the phase diagram should be similar 
to Fig. 2. Furthermore, for the double layer case with spin S = 1/2, we have shown that in 
the net spin model, the critical value is J3 = J4 = J 2 = 0.4288Ji. While other studies of the 
"standard" double layer antiferromagnetic Heisenberg model (corresponding to J3 = J4 = 
in the net spin model), obtained the critical value as J2/J1 = 0.3960 ± 0.0003 (T^j. These 
two values are quite close, suggesting that our results are insensitive to the values of J3 and 
J4. Essential physics is determined by the two couplings: J 2 and J\. 

In summary, we have studied a class of quantum spin models, the net spin model, defined 
by the Hamiltonian Eq. ([I]) for any spin S. The model is controlled by three antiferromag- 
netic couplings: J2/J1, J3/J1, J4/J1, and the measure of the dimensionality L. When 
2J 2 = J3 + J4, the completely dimerized state ipD is an eigenstate of the model and it is 
also the ground state if J\ is sufficiently large. Moreover, wave functions such as if) m ,m' an d 
ipn are also eigenstates of the model. So a subset of the Hilbert space of the net spin model 
is exactly known. To fully understand the model, we have analytically studied the case of 
J2 — J3 — J4 extensively and obtained phase diagrams for both ID (ladder, L = 1) and 
2D(double layer, L = M) cases. For the ladder case, when S — 1/2, L = 1, J 2 = 0.5Ji, 
J3 = Ji, and J4 = the model goes back to the Majumdar- Ghosh model [1]]. Thus, the net 
spin model is a generalization of the MG model for any spin S in two dimensions. We also 
showed qualitative different behavior between ID and 2D, and discussed the crossover from 
ID to 2D. 
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